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Phi meson production in heavy-ion collisions at SIS/GSI energies (~ 2 

'> 

psj ' GeV/nucleon) is studied in the relativistic transport model. We include 

o ' 

contributions from baryon-baryon, pion-baryon, and kaon-antikaon collisions. 

o 

On 



The cross sections for the first two processes are obtained in an one-boson- 
exchange model, while that for the last process is taken to be of Breit-Wigner 
form through the phi meson resonance. The dominant contribution to phi 
meson production in heavy ion collisions at these energies is found to come 
from secondary pion-nucleon collisions. Effects due to medium modifications 

X" 

<-h ' of kaon masses are also studied and are found to reduce the phi meson yield 

by about a factor of two, mainly because of increased phi decay width as 
a result of dropping kaon-antikaon masses. In this case, the <p/K~ ratio is 
about 4%, which is a factor of 2-3 below preliminary experimental data from 
the FOPI collaboration at GSI. Including also the reduction of phi meson 
mass in medium increases this ratio to about 8%, which is then in reasonable 
agreement with the data. 
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I. INTRODUCTION 



There are increasing interests in phi meson production from heavy ion collisions since the 
suggestion by Shor P| that an enhanced production may signal the formation of a quark- 
gluon plasma in these collisions as its production in hadronic matter is usually suppressed 
by the OZI rule [0. However, recent theoretical studies have suggested that because of 
chiral symmetry restoration hadron properties are expected to change in hot dense hadronic 
matter (For recent reviews see Refs. For phi meson in dense matter, QCD 

sum rule calculations have shown that its mass decreases by about 2-3% at normal nuclear 
matter density 0]. Reduced in-medium phi meson mass has also been predicted in vec- 
tor dominance model with dropping kaon masses ||, in hadronic models that include the 
particle-antiparticle vacuum polarization []10f| , and in hidden gauge theory |IT| . A reduced 
phi meson mass in medium is expected to enhance its yield from heavy ion collisions even 
in the absence of the quark-gluon plasma [|12]. In this case, one can also study the phi 
meson in-medium properties via the dilepton spectrum from these collisions as pointed out 



in Refs. [|,|T^,|Tj| . At high temperature, the reduction of phi meson mass may become more 
appreciable because of the decrease of strangeness condensate (ss) as a result of the presence 
of strange hadrons ||. It has thus been suggested that a low mass phi peak besides the 
normal one may appear in the dilepton spectrum from ultrarelativistic heavy ion collisions 
if a first-order phase transition or a slow cross-over between the quark-gluon plasma and the 
hadronic matter occurs in the collisions Jl5| . 

Experimentally, phi meson production has been studied at SPS/CERN by the NA38 



collaboration [T6fl and the HELIOS-3 collaboration |17| using a proton or sulfur projectile at 
200 AGeV and a tungsten or uranium target. A factor of 2 to 3 enhancement in the double 
ratio {4>/{oJ + p°))su{W)l {4 > I{ UJ + P°))pW has been observed. Various theoretical attempts 



have been made to understand this enhancement |12|,|T8|,|T9[ . In particular, an enhancement 
of phi meson yield may be a signature of the formation of a quark gluon plasma in the 
collisions pf. However, the enhancement can also be explained in hadronic models if one 
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takes into account the reduced phi meson mass in medium |12| . In Ref. ||18|| , the color rope 
fragmentation and hadronic rescattering have been introduced to understand the observed 
phi meson enhancement. 

Phi mesons have also been measured at AGS/BNL by the E-802 collaboration in central 
collisions between a 14.6 AGeV/c Si beam and a Au target 2l| . They are identified from 



the invariant mass spectrum of K + K~ pairs, and the measured phi meson mass and width 
are consistent with that in free space. The ratio of the phi meson yield to the K~ yield has 
been found to be about 10%, which can be understood if thermal and chemical equilibrium 



are assumed at freeze out with a temperature of about 110 MeV [21, 22 1 . On the other hand, 



calculations by Baltz and Dover |23j based on the coalescence model, in which the phi meson 
is formed from the kaon and antikaon at freeze out, underestimate the data by a large factor. 
Phi mesons from heavy- ion collisions at SIS / GSI energies are being studied by the FOPI 



collaboration p4| also through the K + K~ invariant mass distribution. About 30 phi mesons 
have been identified in central collisions of Ni+Ni at 1.93 AGeV. Based on these preliminary 
results it has been concluded that the phi meson yield is about 10% of the K~ yield, which 
is very similar to that observed at the AGS energies. This is somewhat surprising since the 
SIS energies are below the phi meson production threshold in the nucleon-nucleon collision 
in free space, while the AGS energies are well above the threshold. 

In this paper, we shall carry out a detailed transport model study on phi meson pro- 
duction in heavy-ion collisions at SIS energies. There are mainly two reasons to work at 
these energies. First, particle production at subthreshold energies is sensitive to medium 
modifications of hadron properties 0J|,^,^]. Secondly, the reaction dynamics of heavy- ion 
collisions at SIS energies is relatively simple compared to that at higher energies as it in- 
volves mainly the nucleon, delta, and pion. In a previous study of phi meson production 
at SIS energies JH], only the contribution from the reaction KK — > cf) has been included, 
and the (f)/K~ ratio is predicted to be less than 1%, which is an order of magnitude smaller 
than the experimental result from the FOPI collaboration. To explain this large discrepancy 
we include here phi meson production also from baryon-baryon and pion-baryon collisions. 



The cross sections for these processes are calculated using an one-boson-exchange model 
with parameters fitted to available experimental data. Some results from this study have 
been reported earlier [p7| . We find that in heavy ion collisions the dominant contribution 
to phi meson production is from the collision, and its inclusion in the transport mode 
brings the predicted <p/K~ ratio closer to the measured one by the FOPI collaboration. 
Furthermore, we study the effect of medium modifications of kaon properties on phi meson 
production. Kaplan and Nelson have pointed out that the kaon mass changes in nuclear 
matter f28j. In chiral perturbation theory, the kaon mass is predicted to increase slightly 
while the antikaon mass decreases significantly as the density increases [f29fl . Overall, the 
total mass of a kaon-antikaon pair is lowered in dense matter. This affects appreciably the 
phi decay width T^kRi which is only about 3.7 MeV in free space, since the phi meson 
mass is very close to twice the kaon mass. A reduced kaon in-medium mass thus increases 
^4>~*kRi so phi mesons are more likely to decay before freeze out, leading to a lower final 
phi meson yield. This reduction is, however, compensated by effects from the dropping phi 
meson in-medium mass. 

The paper is organized as follows: In section [TT] we calculate the elementary cross sections 
for phi meson production and absorption using an one-boson-exchange model. Section fTJ 
describes the nucleon, pion and kaon dynamics in heavy ion collisions using the relativistic 
transport model. Results on the phi meson yield and distributions from these collisions are 
given in section IIV]. Section M is the summary and outlook. 



II. ELEMENTARY CROSS SECTIONS 

For heavy ion collisions at SIS energies, the most abundant particles are the nucleon, A 
resonance, and pion. Phi meson can thus be produced from reactions such as NN — >NN(f), 
NA — >NN<p, AA — >NN<ft, ttN — ><f)N, and 7rA — ><pN. The production of phi meson from 
7T7T scattering is expected to be unimportant due to its higher production threshold than 
that for the other reactions, and is thus ignored. Although both kaon and antikaon are only 
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scarcely produced, phi meson production from the reaction KK — ><p turns out nonnegligible 
in the case of using in-medium kaon masses and will be included. 

Unfortunately, the cross sections for these reactions are not well-known empirically. 
There exist some data for the nN channel near threshold ||30|| , but the NN channel has 
been studied only at high energies [51]. Because of the scarcity of experimental data, an 



one-boson exchange model has been introduced in Ref. |27| to extrapolate from the data 
at high energies to make predictions at energies near the production threshold, and also 
to calculate the cross sections for reactions involving the delta resonance as they cannot be 
measured experimentally. In that study, only the lowest order Born amplitude for phi meson 
production are considered, and this should be a reasonable approximation to the full ampli- 
tude because phi meson production has a very small probability in these reactions. In the 
following, we discuss this model and its predictions for both phi production and absorption 
cross sections. 

The interaction Lagrangians needed in the one-boson exchange model are well-known, 
and they are given by 

ArJW = -^Wf*l> ■ d,7T, (1) 

£ p nn = -g P NNi''l fl rip ■ p M - ^^-^a^Tlf) ■ {d^py - dvPp), (2) 
C p na = z^^VtV 17 • iPuPv - d u p,) + h.c, (3) 
C ^ = ^e^yd-p? • tt. (4) 

In the above, ip is the nucleon field with mass m^; ip^ is the Rarita-Schwinger field for the 
spin-3/2 delta resonance with mass and tt, p M , M are meson fields with masses m n , 
nip, m^, respectively. The isospin matrices, Dirac 7-matrices, and Levi-Civita tensor are 
denoted by r, 7 M , and e^ VOL ^ , respectively. The coupling constants and cutoff parameters in 
£ n NN, £ p nn, and CpNA are taken from the Model II of Table B.l in Ref. P2| . From the 
measured width T^ np w 0.6 MeV, the coupling constant ~ 1.04 is determined. 
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A. Phi meson production from pion-baryon collisions 



The Feynman diagrams for the reactions ttN — ><f)N and ttA — ><j)N are shown in Fig. |T|. 
A monopole form factor is introduced at the 7rp0 vertex, i.e., 
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(A w ,) 2 -m 

The cut-off parameter A£ ^ is adjusted to give a reasonable fit to the experimental data from 
the reaction n~p — > ncf) [jSOfl , and a value of 1.2 GeV is obtained. The comparison between 
calculated results and the experimental data as well as the parameterization by Sibirtsev 
33| can be found in Ref. p7| . 



Since the isospin degree of freedom is usually not explicitly treated in transport model, 
we have thus calculated the isospin-averaged cross sections for ttN — xpN and it A — xpN. 
The results are shown in Fig. |2|. Near threshold the cross section for ttN — ><pN is seen to 
be much larger than the cross section for 7rA — ><f)N as a result of the strong tensor coupling 
at the 7rNN vertex. 

To determine the momentum of produced phi meson in the transport model, we also 
need the differential cross sections for these reactions. Since only two particles are in the 
final state, the magnitude of the phi momentum is fixed by kinematics. The differential 
cross section da/dcosO for these reactions at an energy of 0.2 GeV above the threshold is 
shown in Fig. || 

B. phi meson production from baryon-baryon collisions 

The Feynman diagrams for the reactions NN — >NN(p, NA — >NN<p, and AA — >NN(p 
are shown in Fig. |j. Exchange diagrams are also included as the two final nucleons are 
identical. There is some ambiguity in treating the irpcf) vertex in Fig. ||, which involves two 
virtual particles. Such a vertex does not exist in the Bonn potential model, and there is also 
no theoretically well-defined criteria for prescribing its form. The following form has been 
used in Ref. f27j, 
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2X /, xr V 2 _ m 2 



where g p and are four momenta of the virtual p and 7r meson, respectively. The cut- 
off parameter A^. ^ is the same one introduced earlier in calculating the cross section for 
ttN — >(j)N. The other cut-off parameter A^. ^ is introduced to take into account the vir- 
tuality of the pion. The cross section for pp — >pp<p has been measured at several energies 
|3j| , and the one at the lowest energy has been used in Ref. [27| to fix A 77 and a value 



of 0.95 GeV has been obtained. Comparisons of calculated results with both experimental 



data and the results by Sibirtsev based on an one-pion exchange model |33[ have been given 



in Ref. [ 27j . We note that both A^ and A^. ^ are in the order of f GeV, similar to typical 



values for cut-off parameters in both the Bonn potential model and other hadronic models. 

For the reactions NA — >NN<p and AA — >NN<p, the exchanged pion can be on shell, 
so the cross sections are singular in certain kinematical region. This singularity can be 



regulated by the Peierls method in which the energy of the delta resonance is taken 
to be complex to account for its finite lifetime. As a result, the four-momentum of the 
exchanged pion acquires an imaginary part, which thus renders the pion propagator finite 
at the pion pole. 

The isospin averaged cross sections for phi meson production from baryon-baryon colli- 
sions are shown in Fig. |5|. There are significant differences between the cross sections for 
the reactions NN — >NN<p, NA — >NN<p, and AA — >NN<p. Near threshold, all are much 
smaller than the cross section for phi meson production from pion-nucleon collisions. 

Since there are three particles in the final state in baryon-baryon collisions, the momen- 
tum of produced phi meson is given by the triple differential cross sections d 3 a/dp x dp y dp z or 
double differential cross section d 2 a/p 2 dpdQ if there is an azimuthal symmetry. To simplify 
the numerical procedure, we introduce the assumption that the phi meson momentum is 
isotropic in baryon-baryon center-of-mass frame. We then need only the momentum spec- 
trum da /dp, which is shown in Fig. |] for these reactions again at an energy of 0.2 GeV 
above the threshold. 



C. Phi meson production from kaon-antikaon annihilation 



The total cross section for KK — xj) is assumed to have a Breit-Wigner form as in Ref. 



Tl 



i.e., 



3_7r (m^)' 

k 2 (M 2 - m 2 ) 2 + {m+rtf 



a(KK ><f>) = Z7*?, m .r.» > ^ 



where is the phi meson decay width to KK and is given by 

r * " e^ 2 ' (8) 

with the coupling constant gL^-/47r rs 1.69 determined from the empirical width of 3.7 
MeV. In using Eq. (0), we have thus neglected the small effect due to the decay of phi 
meson into other channels. 



D. Phi meson scattering by nucleons 

Once a phi meson is produced, there are three kinds of reactions it can take part in: 
decay, inelastic and elastic scattering by nucleons. The dominant decay mode of phi meson is 
— > KK, with a decay width given by Eq. (^). Here we only consider phi meson absorption 
due to <pN collisions as that due to (pir and <pA scattering is expected to be insignificant 
because of the low pion and delta densities compared with that of nucleons (see Fig. |8|). We 
have included the reactions <j)N — >KA, <pN — >pN, <pN — >tcN, and <j)N — >7rA. The cross 
section for the first reaction is related via detailed balance to that of AK — * <f)N, which has 
already been calculated in Ref. |12[ based on a kaon exchange model with parameters taken 



from Ref. ||35|| . The reaction (fiN — > KT, is neglected since the NUK coupling constant 
(<7jv£K-/47r ~ 0.6) is much smaller than that of NAK (g 2 N \ K /^ ~ 16) [j3~5[j . The cross 
section for the other three reactions can be obtained from those for pN — > (pN, ttN — ► <pN, 
and 7rA — > <pN using detailed balance relations. We have already determined in the above 
the cross sections for phi meson production from pion-baryon collisions. For the reaction 
pN — > <pN, we simply interchange the pion and rho meson in the first diagram of Fig. [l|. 
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The cross sections for these phi meson absorption processes are summarized in Fig. |7]. We 
see that 4>N — >KA is the dominant one as already pointed out in Refs. p|,|P2||, because this 
is the only one among the four reactions that is not suppressed by the OZI rule . 

The phi-nucleon elastic scattering is also included. Experimentally, the cross section for 
this process can be extracted from data on phi meson photoproduction using the vector 
meson dominance model |36|| . In the energy range 3.5 GeV< pi a b < 5.8 GeV, a value of 



0.56 mb has been obtained. Here we make the rough approximation that the cross section 
is a constant of 0.56 mb in the energy range considered. We note that a slightly larger value 



has been obtained in Ref. []37|. The phi-nucleon scattering does not change the phi meson 
yield, but affects its momentum distribution. Since this cross section is rather small, the 
effect is not appreciable as shown below. 

E. Discussions 

So far we have calculated exclusive phi meson production cross sections with no pions in 
the final state. In principle, phi meson can also be produced in association with one or more 
pions, and one needs for heavy-ion collisions inclusive cross sections such as NN — > (ftX. 
Since the beam energies we consider in this work are below the phi meson production 
threshold in nucleon-nucleon collision, we expect that exclusive reactions such as NN — > 
NNtfi, which have lower thresholds, are more important than the ones with pions in the final 
state, which require additional energies and are thus suppressed. 

We would like to point out that in spite of the success of the one-boson-exchange model 
at low energies, this model has not been carefully verified at energies above the phi meson 
production threshold. It is often believed that the quark degree of freedom becomes relevant 
at a kinetic energy of 1-2 GeV, so it is not known if our procedure of using the data 
point at pi a b = 10 GeV to fix the cross section for pp — >pp<ft near the threshold is reliable. 
Nevertheless, this model has been used quite often in calculating the cross sections for 
pion production P5[ , eta production pSfl , kaon production [|I],[|T[], and baryon resonance 
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excitations [42] from nucleon-nucleon collisions in this energy region. Measurements of 



pp — > ppcj) at low energies are being carried out at SATURNE [J3J, and these data will be 
very useful in testing the validity of our model. With this caveat in mind, we proceed to 
our transport model study using the above elementary cross sections. 



III. THE NUCLEON, PION AND KAON DYNAMICS 



A. The relativistic transport model 

For the description of heavy-ion collisions at SIS energies, we use the relativistic transport 
model (RVUU) developed in Ref. [P] from the non-linear o — uj model and extensively 



used in previous studies. At these energies, the colliding system consists mainly of nucleons, 
delta resonances, and pions. While medium effects on pions are neglected as in most trans- 
port models, nucleons and delta resonances propagate in a common mean-field potential 
according to Hamilton equations of motion, i.e., 

where E* = (m* 2 + p* 2 ) 1 / 2 . The effective mass and kinetic momentum of a baryon are given, 
respectively, by m* = m — g a (a) and p* = p — guj(uj), where the expectation values (a) 
and (uj) are related to the attractive scalar potential and the vector current of a baryon 
in nuclear matter. We use in this work the parameter set that corresponds to the soft 



equation of state [^6J with a compressibility K = 200 MeV and an effective nucleon mass 
m* N = 0.83 m^v at normal nuclear density. These particles also undergo stochastic two-body 
collisions, including both elastic (NN NN, NA -> NA, AA -> AA) and inelastic 
(NN <-» NA, A <-> Nn) scattering. The standard Cugnon parameterization |57j and proper 
detailed-balance prescription ||8| are used for describing these reactions. 

Besides pions, other particles are also produced in heavy-ion collisions. However, their 
production probabilities are very small at these energies, and they can thus be treated 
perturbatively, in the sense that their production and interactions do not affect the dynamics 
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of nucleons, delta resonances, and pions. This is carried out using the perturbative test 



particle method of Ref . P9| , |50f . In the case of kaons, they are produced from pion-baryon 
and baryon-baryon collisions whenever the energy is above the threshold. Then, a probability 
factor given by the kaon production cross section to the total pion-baryon or baryon-baryon 
cross section is assigned to the produced kaon. The motion and collision of the kaon are then 
followed. However, only the kaon momentum is changed after a kaon-nucleon collision. This 



method allows us to include easily the medium effects on the kaon as shown in Ref. |49 



We use in this study the kaon in-medium mass obtained from the mean-field approximation 
to chiral Lagrangian ||28|| . As in Ref. [51 ] , the kaon mass m* K at nuclear matter density p is 
given by 

1 I 1 / 2 



m K = m K 



1 



,3 p 



8 J 2 ' 



(10) 



In the above, ttik is the kaon mass in free space; / ~ 93 MeV is the pion decay constant; 
and p s is the nuclear scalar density. The kaon-nucleon sigma term is denoted by T^xn and 
is taken to be about 350 MeV. At normal nuclear density, the kaon mass is then increased 
by about 10 MeV. For antikaon, its vector interaction is opposite to that of kaon, so the last 
term in Eq. ( |T0"P has instead a negative sign. The antikaon mass at normal nuclear density 
is reduced by about 100 MeV. Corrections to the mean-field results may not be negligible 
29| , |52|| and are still under debate. Since there are other uncertainties in the treatment of 



kaons such as their production cross sections from baryon-baryon collisions, we will not 
address the sensitivity of our results on these corrections. 



Our treatment of kaon follows essentially Refs. |50| , [5T| , |53| , where the baryon-baryon pro- 
duction channels BB — > NAK and BB — > NUK are included with their cross sections 



given by the Randrup and Ko parameterization ||54|| . We have also included the pion-baryon 
channels ttN — > AK and ttN — > Y,K. For these cross sections, the parameterizations due to 



Cugnon []55| are used. Kaon production from ttA collisions is neglected as its cross section 



has been shown to be much smaller than that for nN collisions [BB 



The production of antikaon in heavy-ion collisions at SIS energies was first studied in Ref. 



11 



57 1 using the relativistic transport model and the Zwermann-Schiirmann [5l| parameteriza- 



tion for the antikaon production cross section from the nucleon-nucleon collision. However, 
it has been recently found that this parameterization overestimates the cross sections near 
the threshold by about one order of magnitude ]5^,f)(J 5 an d contributions from pion-hyperon 
interactions [BTI], which were neglected in Ref. [|57|, are needed to account for the measured 



K yield. In this work, we have thus modified our earlier treatment of antikaon produc- 



tion by using instead the parameterizations proposed in Ref. |59j for antikaon production 
cross sections in both pion-baryon and baryon-baryon collisions. We have also included 
contributions from the pion-hyperon collisions. 

B. The pion and nucleon yield and distributions 

In the upper window of Fig. [8], we show the time evolution of central baryon and pion 
densities in Ni+Ni collisions at an energy of 1.93 AGeV and zero impact parameter. It is 
seen that the maximum baryon density reached in the central region is around three times 
normal nuclear density, while the maximum delta and pion densities are about 1/6 of the 
maximum nucleon density. The abundance of pions and deltas is given in the lower window 
of Fig. which shows that the maximum number of delta is around 18. Since the number 
of pions reaches around 30 at freeze out, the pion/nucleon ratio is thus about 0.25, which is 
in agreement with the experimental data [^] as well as other transport model results ||63|| . 

Our results for proton and 7r~ rapidity distributions are compared with the experimental 
data from the FOPI collaboration 0] in Fig. |9|, where y = y / y P rojectHe is the normalized 
center-of-mass rapidity. The theoretical results are obtained for the impact parameter range 
b < 2 fm in order to compare appropriately with the data which have a centrality selection 
of 100 mb. Also, only protons that have a local density less than 0.02 fm -3 at freeze out 
are included. Otherwise, they are considered as bounded in the deuteron and other light 
fragments. The agreement between theoretical results and data, especially in mid-rapidity 
region, is quite good, so we feel that our transport model describes properly the dynamics 
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of nucleons, deltas, and pions. The proton rapidity distribution shows a peak at y = 0, thus 
indicating that considerable amount of stopping is achieved in the collisions. 

In Fig. [10], the nucleon and pion transverse mass distributions at mid-rapidity region 
(—0.2 < y < 0.2) are shown. We find that the inverse slope parameter for nucleons is about 
180 MeV, while that of pions is about 120 MeV. The difference is due to radial flow, which 
has a larger effect on nucleons than pions. 



C. The kaon yield and distribution 



The abundance of kaons and antikaons is shown in Fig. [LI]. Without kaon medium 
effects the K + number reaches about 0.18, while that of K~ reaches about 1.4 x 10~ 3 . 
Including kaon medium effects, the kaon yield is reduced by about 20% to about 0.15, while 
the K~ yield increases by about a factor of 3 to about 4 x 10~ 3 . In our model, K + feels 
a weak repulsive potential in nuclear matter, which increases its production threshold and 
thus reduces the yield slightly. On the other hand, K~ has a strong attractive potential, so 
its production threshold is reduced and a much larger yield is obtained. 

The kaon and antikaon rapidity distributions for central Ni+Ni collisions are shown in 



Fig. |12| for both cases with and without kaon medium effects. As in the case of pions 



and nucleons, the rapidity distributions peak at y — 0. In Fig. [L3|, the kaon and antikaon 
transverse mass distributions in mid-rapidity region (—0.2 < y < 0.2) are shown. For kaon, 
the inverse slope parameter is found to be about 100 MeV without kaon medium effects, and 
increases to about 120 MeV after kaon medium effects are included. The kaon momentum 
thus increases as a result of propagation in the repulsive mean-field potential. For antikaon, 
an inverse slope parameter of about 110 MeV is found in the case without kaon medium 
effects. This is somewhat larger than that of the kaon transverse mass spectra, mainly due 
to the absorption of low-momentum antikaons by nucleons because of a larger absorption 
cross section [[ST]]. Including kaon medium effects, the inverse slope parameter decreases 
to about 100 MeV. The propagation in the mean-field potential thus reduces the antikaon 



13 



momentum as a result of the attractive potential. Therefore, from the ratio of kaon and 
antikaon transverse mass spectra it is possible to learn about the kaon medium effects as 
pointed out in Ref. [[>4[ . 



IV. PHI MESON PRODUCTION FROM HEAVY ION COLLISIONS 



A. Transport model for the phi meson 



Phi mesons are also treated in our model by the perturbative test particle method |49 
Whenever a pion-baryon or a baryon-baryon collision is above the threshold for phi meson 
production, a phi meson is produced and is assigned a probability factor that is given by the 
phi meson production cross section to the total pion-baryon or baryon-baryon cross section. 
The motion of the phi meson, including its collisions with other nucleons, is then followed. 
As in the case of kaon, effects of these collisions are included for phi mesons but not nucleons. 
The numerical simulation of heavy ion collisions is stopped when the number of collisions 
is small. The number of phi mesons at this time is taken to be the yield accessible to 
measurement by identifying KK pairs. As shown in Refs. [50,57], KN and KN scattering 



are significant before freeze out, so the A' A' pair coming from the phi meson decay inside the 
fireball loose their correlation before they get out and cannot be identified as a phi meson 
from the KK invariant mass spectrum. 

In 7iN or 7rA channels phi meson is produced in a two-body final state, so the magnitude 
of its momentum is fixed by kinematics. The polar angle of the phi meson momentum 
depends on the differential cross section and is thus determined by a probability function 
given by the ratio of the differential cross section to the total cross section. For the azimuthal 
angle of the momentum, it has a uniform distribution between zero and 2n. 

The baryon-baryon channels are more complex as the phi meson is produced in a three- 
body final state. In this case, the final momentum of the phi meson follows a multivariable 
probability function given by the ratio of the differential cross section to the total cross 
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section. To simplify the numerical procedure, we assume that the phi meson is produced 
isotropically in the baryon-baryon center of mass frame as discussed in Section [11 B| . This 
assumption is often used in transport model study when no data exist for the differential 
cross section. The magnitude of phi meson momentum is then determined by the momentum 
spectrum such as that shown in Fig. ^] 

When including medium effects, we only change the available phase space due to the 
modification of hadron masses. Modifications on coupling constants, cut-off parameters, and 
the functional form of cross sections are neglected in this study. In the case of a reaction 
A + B — >• C + D, the functional form of the cross section as a function of E(= \/s* — v^o) 
is assumed to be the same as in free space, where the modified center of mass y/s* energy 



is given by y/s* = \Jp*2 + + \]Pb + m B> while the modified production threshold is 
y/s% = m*Q + Trip , with p* and m* being the momentum and mass of a particle in nuclear 
medium. 

In our study medium modifications of pion properties are neglected, while those on 
nucleons are always included. For the treatment of kaon, antikaon and phi meson production, 
we will consider three scenarios. In the first scenario, we neglect medium effects on both 
kaons and phi mesons. It is already known from previous studies of antikaon production 



57| . |60|| that the neglect of kaon medium effects underestimates the antikaon yield. In the 
second scenario, we include medium effects on kaons but neglect those on phi mesons. 
The kaon medium effects are obtained from the chiral perturbation calculation as discussed 
earlier. Inclusion of kaon medium effects improves the agreement of the antikaon yield with 
the experimental data as discussed in Section [111 C|. Since a phi meson decays mainly into a 



kaon-antikaon pair, medium modifications of kaon properties affect the final phi meson yield 
determined from the reconstruction of K + K~ pairs. In the last scenario, we include both 
kaon and phi meson medium effects. For the in-medium phi meson mass we use the results 
from QCD sum rule calculations by Hatsuda and Lee ||], 

« 1.0 - 0.0225— , 11 
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which is obtained by taking the nucleon strangeness content to be about 0.15. 



B. The phi meson yield 

The results for phi meson yield in central Ni+Ni collisions at 1.93 AGeV are shown in 



Fig. [14]. Phi mesons are mostly produced during the first 10 fm/c, so the yield increases 
with time. Afterwards, phi meson decay and reabsorption become more important, and 
the phi meson yield is seen to decrease with time. We find that at t ~ 22 fm/c two-body 
scattering becomes rather scarce, so we take this as the freeze-out time and determine the 
final phi meson yield by the number of phi mesons which have not decayed. 

In the scenario that both kaon and phi meson medium effects are neglected, the phi 
meson yield has a value of 3.3 x 10 -4 , while that of K~ is about 1.4 x 10 -3 . The <p/K~ 
ratio in central collision is thus about 25% and is about a factor of 2 too large compared 



with preliminary FOPI data. However, as pointed out earlier |57],|6(J, the K~ yield without 
medium effects is inconsistent with experimental results. In order to explain the observed 
K~ yield, kaon medium effects are needed, and this will increase K~ yield by about a factor 
of 3. The <p/K~ ratio would then be about 8% and is in reasonable agreement with the 
FOPI data. 

On the other hand, the inclusion of kaon medium effects inevitably affects the phi decay 
width. With a lower in-medium KK threshold, the phi meson decay width increases signifi- 
cantly |14| i so the probability for a phi meson to decay inside the dense matter is enhanced. 



Because of strong final-state interactions of kaons and antikaons, these phi mesons cannot 
be identified from K + K~ reconstruction. Including this effect we find that the final phi 
meson yield is reduced by about a factor of 2 to about 1.6 x 10~ 4 . The <f)/K~ ratio is now 
about 4% and is about a factor of three too small compared to the FOPI data. 

Including also the dropping of phi meson mass in dense matter leads to an increase 
of phi meson yield, as a result of lower production threshold. For example, around 2po 
where most phi mesons are produced, the reduction in the production threshold is about 45 
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MeV. The final phi meson yield after including both kaon and phi meson medium effects is 
about 3 x 10~ 4 , so the (f)/K~ ratio is brought back to about 8%, and is again in reasonable 
agreement with preliminary data from the FOPI collaboration. 

In Figs. |1|, [L6|, and [IT], the phi meson yield from different production channels is shown 
for the three scenarios mentioned above. In all three cases the irN channel is the dominant 
one, contributing about 60% of the total phi meson yield, The importance of the irN channel 
is due to its large cross section near the threshold, although pions have a lower abundance 
than that of nucleons and the production threshold for this process is also the highest. The 
baryon-baryon channels add up to about 15-20% of the total phi meson yield, with NA 
contributing the most. If phi meson production cross sections in NA and A A collisions are 
assumed to be the same as that in NN collisions, we find that the final phi meson yield is 
reduced by about 10%, with the baryon-baryon contribution being only about 5-10% of the 
total. In the scenario that both kaon and phi meson medium effects are neglected, the KK 
channel contributes only about 7% of the total phi meson yield. Inclusion of kaon medium 
effects increases it slightly to about 12% of total phi meson yield. 

The impact parameter dependence of phi meson yield in Ni+Ni collisions at 1.93 AGeV 
is shown in Fig. [Tj| From this, the minimum-biased phi meson production cross section can 
be estimated. It is about 0.15 mb when both kaon and phi medium effects are neglected. 
This reduces to about 0.05 mb if kaon medium effects are included but those on phi mesons 
are neglected. Including both kaon and phi meson medium effects increases the production 
cross section to about 0.13 mb. 

Both phi meson decay and absorption have significant effects on its final yield. Without 
kaon medium effects, about 20% of produced phi mesons are absorbed in cf)N scattering, 
and about another 20% decay before freeze out. Including kaon medium effects, absorp- 
tion by nucleons still removes about 20% of produced phi mesons, but phi meson decay is 
significantly enhanced, leading to more than 50% reduction of the phi meson yield. 
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C. The phi meson rapidity distribution and transverse mass spectrum 



The rapidity distribution and transverse mass spectrum of phi mesons at freeze out are 



shown in Fig. 19 and Fig. g(J respectively. The phi meson rapidity distribution is seen to 
peak at y = as those of nucleons, pions, and kaons. Their transverse mass distribution has 
a typical exponential form with an inverse slope parameter of about 110 MeV for all three 
scenarios considered. 

The results shown in the above include (f)N scattering with a cross section of 0.56 mb. 
To see more clearly the effects due to scattering we show in Fig. [H] the phi meson transverse 
mass spectra with and without <pN elastic scattering in the scenario of no medium effects 
on kaons and phi mesons. To facilitate the comparison, different results are normalized to 
have the same value at m t — = 25 MeV. We find that the inverse slope parameter for phi 
mesons in the mid-rapidity region, obtained by fitting the spectrum up to 0.4 GeV, changes 
only slightly when rescattering is included. As the extraction of <pN elastic cross section 
from the experimental data is model-dependent, its value may not be well determined. To 
see the effects due to a different cross section, we have also repeated the calculation using 
a <pN elastic cross section of 8.3 mb, which is the total (f)N cross section estimated from 



the phi meson photoproduction data |65j and serves as an upper bound on the elastic cross 
section. The inverse slope parameter in this calculation is found to be about 130 MeV. 
The increase of the slope parameter due to scattering can be understood as follows. At SIS 
energies most produced phi mesons have low momenta, so scattering with nucleons increases 
their momenta and helps them in achieving thermalization with nucleons. Similar effects 
have been found for kaons in heavy ion collisions at subthreshold energies [|50| ,|66|1 . 

V. SUMMARY AND OUTLOOK 

We have carried out a transport model study of phi meson production from heavy ion 
collisions at SIS energies. The production channels included in our study are ttN — >4>N, 
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ttA — xf)N, NN — >NN(j>, NA — >NN<j), AA — >NN(j), and KK — ><j). The cross sections 
for the first five reactions are obtained from the one-boson-exchange model, while that for 
the last one is taken to be of Breit-Wigner form through the phi meson resonance. More 
than half of the phi mesons are found to be produced from 7rN collisions. 

The predicted 4>/K~ ratio depends on whether medium effects are included for kaons and 
phi mesons. Neglecting both kaon and phi medium effects, the ratio is about 25%, which 
is about a factor of two larger than that from the FOPI data. However, the K~ yield in 
this case is about a factor of 3 below the measured value. To account for the observed K~ 
yield requires the inclusion of medium effects on kaons, such as that predicted by the chiral 
perturbation theory, so that the antikaon production threshold is reduced. The 4>/K~ ratio 
turns out to be about a factor of 2-3 smaller than the experimental results. Including also 
dropping phi meson mass in nuclear medium as predicted by the QCD sum rules increases 
the phi meson yield. The final <j)/K~ ratio becomes about 8% and agrees reasonably with 
the preliminary FOPI data. 

It will be of interest to also measure phi mesons through their dilepton decay channel, 
which will be carried out in the near future by the HADES collaboration Jj7[ at GSI. 
Since both the KK and the dilepton channels are subject to same uncertainties in the 
production mechanism, the ratio of the phi meson yield measured from the KK channel to 
that measured from the dilepton channel is expected to be less sensitive to such uncertainties. 
In addition, the dilepton channel can also reveal directly the properties of vector mesons 
in nuclear medium, especially the change in their masses |63|j68|j69[| . For phi mesons, this 



is not possible via the KK channel since the correlation of a KK pair is lost due to KN 
scattering in the nuclear medium. Furthermore, the shape of the phi meson peak in the 
dilepton spectrum can also indirectly give information on the kaon medium effects |70| . For 
this purpose, a careful study of the background is needed to determine whether the phi 
meson peak will still be observable. Previous study |63| shows that the background due 



to nucleon and pion bremsstrahlung as well as the eta Dalitz decay is small near the phi 
meson peak. The major background is thus expected to come from decays of rho and omega 
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mesons. A careful study of the production of vector mesons in heavy ion collisions is needed 
in order to learn the phi meson properties from the dilepton decay channel. Such a study is 
in progress. 
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FIGURES 



n ^n.^^-^^ <P 

FIG 1. Feynman diagrams for ttN — >(f>N and irA — ><j>N. 
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FIG. 2. Isospin-averaged cross sections for 7r7V — >4>N (solid curve) and 7rA — >4>N (dotted 
curve) . 
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FIG. 3. Differential cross sections for ttN — ><fiN (solid curve) and 7rA — ><fiN (dashed curve) 
at 0.2 GeV above the threshold. 
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FIG. 4. Feynman diagrams for NN — >NN(j), NA — >NN(p, and AA — >NNcp. 
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FIG. 5. Isospin-averaged cross sections for NN — >NN(p (solid curve), NA — >NN(f> (dashed 
curve), and A A — >NNcf) (dotted curve). 
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FIG. 6. Phi meson momentum spectrum from the reactions NN — >NN<f>, NA — >NN<f>, and 
AA — >NN4> at 0.2 GeV above the threshold. 
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FIG. 7. Cross sections for 4>N — >Kk (solid curve), 4>N — >7riV(dotted curve), 4>N — >pN 
(dashed curve), and 0iV — >ttA (dash-dotted curve). 
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FIG. 8. Upper window: the time evolution of central densities of baryons (solid curve), deltas 
(dashed curve), and pions (dotted curve) in Ni+Ni collisions at 1.93 AGeV and impact parameter 
fm. Lower window: the time evolution of the delta (dashed curve) and pion (dotted curve) 
abundance. 
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FIG. 9. The rapidity distribution of nucleons and pions in Ni+Ni collisions at 1.93 AGeV and 
impact parameter b < 2 fm. 
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FIG. 10. Same as Fig. |9| for the transverse mass distribution. 
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FIG. 11. The abundance of K + and K~ with (solid curve) and without (dotted curve) medium 
effects in the same reaction as in Fig. |8|. 
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FIG. 12. The rapidity distribution of K+ and K~ in central Ni+Ni collisions at 1.93 AGeV 
with (solid curve) and without kaon (dotted curve) medium effects. 
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FIG. 13. Same as Fig. [l2| for the transverse mass distribution of K + and K~ with (solid curve) 
and without (dotted curve) medium effects. 
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FIG. 14. The phi meson yield in Ni+Ni collision at 1.93 AGeV and impact parameter of fm 
for three different scenarios. 
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FIG. 15. The time evolution of phi meson yield from different channels in central Ni+Ni colli- 
sions at 1.93 AGeV without kaon and phi meson medium effects. 
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FIG. 16. Same as Fig. 
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with kaon medium effects but no phi meson medium effects. 
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FIG. 17. Same as Fig. 
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with both kaon and phi meson medium effects. 
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FIG. 18. The phi meson yield in Ni+Ni collision at 1.93 AGeV at different impact parameters 
for three different scenarios. 




FIG. 19. The rapidity distributions of phi mesons in central Ni+Ni collisions at 1.93 AGeV for 
three different scenarios. 
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FIG. 20. Same as Fig. ^ for the transverse mass distribution. 
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FIG. 21. The transverse mass distribution of phi mesons in central Ni+Ni collisions at 1.93 
AGeV with (dotted curve for cr^v = 0-56 mb and solid curve for cr^v = 8.3 mb) and without 
(dashed curve) <f>N scattering. 



